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Abstract 

A multiple access relay channel (MARC) is considered in which an analogue-like network coding is implemented 
in the relay node. This analogue coding is a simple addition of the received signals at the relay node. Using "nulling 
£f) \ detection" structure employed in V-BLAST receiver, we propose a detection scheme in the destination which is able 

to provide a diversity order of two for all users. We analytically evaluate the performance of our proposed scheme 
for the MARC with two users where tight upper bounds for both uncoded and Convolutionally coded transmission 
blocks are provided. We verify our analytical evaluations by simulations and compare the results with those of non- 
cooperative transmission and Alamouti's scheme for the same power and rate transmission. Our results indicate that 
while our proposed scheme shows a comparable performance compared to the Alamouti's scheme, it substantially 
outperforms the non-cooperate transmission. 

Index Terms 

Cooperative communication, MARC, analogue network coding, V-BLAST detection, Alamouti space-time coding 

I. Introduction 

It has been shown that a cooperative communication can provide more reliable communication by reducing the 
00 ■ error probability, while the required power remains the same. The concept of cooperation in the communication 
^ , has a long history going back to £0 where an information theoretic aspect view of cooperative communication in 
q \ a multiple access channel (MAC) is defined. Although this work was based on the concept of feedback, it can be 
^ ■ seen as a beginning of the notion of cooperation in communication. This concept have then been further extended 
O , in the information theoretic framework in several seminal works in (21-||71. 

The initial spark started in |@] then lead to a fundamental context of network coding (NC), as a method for 
T"! • achieving capacity of wired networks, introduced in 2000 [HI. Thereafter, there have been a vast number of researches 
.£h ' on trus ^ ea w here different advantages of using NC have been addressed. Owing to the cooperative nature of 
^ | network coding, in [9], it was shown that using NC over a wireless network can provide diversity, where it is 
^ ■ called cooperative diversity in NC, similar to its previous counterparts iTTOl — ITT31 . It then lead the researchers in 
wireless communication also to become interested in NC, either in ad-hoc or cellular networks ||T6l - liT9"1 . Moreover, 
implementing NC, as an alternative method of routing, in which the nodes in the network process their received 
packets, can help to increase the multiplexing gain of such networks. 

A very simple example of a cooperative system is a multiple access relay channel (MARC). Although there is 
no direct cooperation between the users, sharing the relay node can be seen as a cooperative behavior. A relay 
node can be utilized conventionally, where at each time slot the relay only serves for one of the users. On the 
other hand, if the relay can process the received data from the users and aggregate these data and then send the 
new aggregated data to the destination, it would be possible to achieve diversity. One instance of such doing can 
be found in EOl . where an XORing scheme was considered as a simple way of implementing NC in a wireless 
network in order to provide diversity in the network. The authors in [21] have proposed a scheme in a MARC in 
which the users send the log-likelihood ratio (LLR) of their coded messages to the destination in a noisy channel, 
and they have shown that their proposed scheme improves the bit error probability (BEP) of the system. Practical 
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analogue NC for cooperation in the physical-layer of MARC has not been considered yet, and it has prompted us 
to consider this method of cooperative scheme here. 

In this paper, we propose a cooperative scheme based on a simple addition of the received signals in the relay 
node. We loosely call this scheme as an instance of analogue NC. To this end, we consider a system which uses 
Decode and Forward (DF) for relay cooperation, and after modulating the users' re-encoded packets, simply adds 
the packets for transmission in the relay. We assume that each node has a single antenna. Also we assume that 
each user has a power limitation of P , but the power limitation of the relay node in n-user MARC will be nP 
which is logical. We propose a detection method in the end-node, i.e. destination, being based on V-BLAST model 
in a multiple-antenna system 11221 . fl23l . We derive a tight upper bound on the system's BEP for both uncoded 
and coded transmissions, from which we show that using V-BLAST detection technique, we can achieve full-order 
diversity, being two in a two-user MARC setting. Simulation results indicate that the bounds are tight enough. 
Moreover, results show that the proposed scheme can achieve a diversity order of two for any number of users and 
only one relay node , in contrast to the traditional methods (e.g. Alamouti scheme) in which each user must have 
two antennas for transmission. 

The paper is organized as follows. In Section JD we briefly describe the system and the proposed transmission 
protocol. In Section [III] we introduce the detection algorithm. The performance analysis is presented in Section HVl 
Numerical experiments are shown in Section E3 We finally finish the paper with the concluding remarks in 
Section ED 

Notation and abbreviations in the paper are as below. Boldface lower case letters denote column vectors. 
Concatenation of several vectors is also denoted by a boldface lower case letter, k — th element of the vector 
x is denoted by x^. Boldface upper case letters are used to denote matrices. . T denotes transpose operators. The 
cadinality of a set S is shown by \S\. Pr(A) denotes the probability of event A. When dealing with random 
variable, x, we use Fx{x), fx(x), and <p x (S) to denote the cumulative distribution function (CDF), probability 
density function (PDF), and characteristic function, respectively. Finally, the notation E x [gi(a;)] is used to denote 
taking expectation of g(x) with respect to the subscript x. 

II. System Description 

A. System Model 

For the simplicity of presentation, we consider MARC with two users A and B, and one relay R, and the 
destination D. The results can be generalized for more than two users sharing one relay. Figure [TJ shows a two-user 
MARC used throughly in this paper. 




Fig. 1. An illustrative view of a multiple access relay channel with two users. 

In what follows, we model the received signals at the destination and relay node. The transmitted signal by the 
relay node is considered in the next subsection. Let /i t r denote the channel gain between nodes t and r, where 
t G {A,B,R} and r £ {R,D}. The channel gains are assumed to be independent zero mean complex Gaussian 
random variables with the powers of SI^r and S1_>d> for the link between the user (A or B) and the relay (R), 
and the link between node (A, B, or R) and the destination, respectively. That is, we have considered different 
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qualities for user-relay and user (relay)-destination channels. We assume a block fading in which the channel (gain) 
coefficient is constant during each block transmission ll24l . Moreover, we assume that each user sends a codeword 
Xj G X. 

For all the wireless communication links shown in Figure Q] we consider an additive white Gaussian noise 
(AWGN) with one-sided power spectral density (PSD) of Nq. Then, knowing the fact that hf^'s are normally 
distributed with the powers mentioned above, the instantaneous signal to noise ratio (SNR) of the link between 
user t G {A,B} and the relay node R is given by SNR^ = jj-P\ht,ii\ 2 being distributed exponentially as follows 



SNR tR ~exp(--N n^ R ) 
Furthermore, the received signals at the destination from two users can be modeled as 



Ybd 



Vdx a + n A 

Vd-Xb + n E 



(1) 

(2) 
(3) 



whereas ni and ri2 are the received additive (vector-valued) white Gaussian noises at the destination side. The 
x B and x A , the vectors of dimension n, denote the transmitted blocks from users A and B, respectively, with the 
following power constraints: 



1 - 



P 



P 



i=l 



where n is the number of bits in each block. Similarly, the received signals at the relay node are given byS 

Yar = /ia,rx a + n A 
Ybr = Vr-Xb + n B 

where n A and rig are the noises received at the relay node through user A and B, respectively. 



(4) 
(5) 



(6) 
(7) 



B. Transmission Algorithm 

Since, in this paper we aim at proposing a cooperative protocol while sharing a single relay node by two nodes, 
we propose to divide the time as shown in Figure [2] In Figure [2] the transmission time slots are allocated based on 
time division multiple access (TDMA) scheme. As illustrated in Figure |2j in the first two time slots, users A and 
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Fig. 2. The transmission time slots allocated for users and relay node. 



B transmit their blocks. Then, the relay, based on its received signals y A R and y B R, decodes the received blocks. 
Depending on whether or not the relay can correctly decode the received blocks, which can be recognized by a 
error-detecting code (e.g., CRC codes), the relay decides upon cooperation with the two transmitters in the third 
slot, as shown in Figure |2] leading to four possible scenarios summarized in Table U 

As can be inferred from Table [TJ relay cooperation can be categorized into four states: First, corresponding to So 
, the relay will be silent, and as a result the destination only receives the direct block transmissions from two users. 
Hence, the bit error probability is the same as that of a direct transmission with power P. Second, corresponding 
to si and S2, the relay node only cooperates in transmission of one of the users. Thereby, only one of the users 

'To keep generality, we assume different statistics for the additive noises corresponding to users' signals received at the relay and destination 
nodes. 
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Event 


State Notation 


Relay Action 


Relay cannot decode none of its received blocks 


so 


Silence 


Relay can decode only received block of user A 


Sl 


Transmitting xa 


Relay can decode only received block of user B 


S2 


Transmitting xb 


Relay can decode the received blocks of both users 


S3 


Transmitting xa + xb 



TABLE I 

Four possible states at the relay node. 



benefits a second order diversity transmission. Denoting this user by u (u £ {A,B}), the received signal at the 
destination from the relay will be 



(8) 



where x u and are relay's estimate of user's signal and additive noise at the destination, respectively. In this 
case, the bit error probability for the assisted user is equal to that of a Maximum Ratio Combining (MRC) receiver, 
11251 , with two independent receptions of the transmitted signal, and for the other user, it is similar to the first state. 
Lastly, for S3 , the relay can successfully decode both users' blocks. Then, the relay encodes and modulates these 
signals separately, and transmits the sum of the modulated analog signals. Thus, the destination node will receive 
the following signal 



Yrd = Vd(x a + x B ) + n R . 



(9) 



III. V-BLAST-Based Detection Algorithm 

In this section, we propose a method for detection at the destination node. First, we assume that the destination 
knows fading coefficients, e.g. by acquiring transmitting pilot symbols, and also knows the state of cooperation at 
the relay, e.g. it can be implemented by sending 2 bits encoding four states. For so in Table H both users have 
interference-free direct transmission, and, as stated above, the conventional receiver is used. For si (S2), the user 
which is able to utilize relay resource, the destination uses MRC of the two received signals, one from the user 
and the other from the relay. In the last state, we propose a novel method for detection relying on the basis of 
V-BLAST detection. First, from equations we rewrite the three received signals in three time slots (see 

Figure [2]) at the destination in a following matrix format given by 




h 



h 



A,D ) 




f 







^-B,D 


R,D J 







x b + 



(10) 



Without loss of generality, we focus on detecting user A's block. Doing so, we first define the vector conveying 
user B's as follows 



^B,D 2 + ^R,D 2 




(11) 



The null space of u B , denoted by U = [ui 112], can then be spanned by 





u 2 



-h* 



h 



B,D 




(12) 



It should be noted that, the matrix U defined in equation (112] ) is one of many possible null spaces for vector 
Ua- Now multiplying both sides of (fTOt by transposed of matrix U, one may obtain 





( Yad ^ 






Ybd 






V Yrd ) 





/iR.D^B.D 



y/\h B , D \ 2 + \h R . D \ 2 
^A,D 



X A + 



=n 



flB.D 



VI^b.d| 2 +|/ir.d| 2 B ^\h B , D \ 2 + \hR,D\ 



-n'i 



< 



(13) 
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As can be realized from equation ([T3T i, we have two receptions for user A's signal. Owing to independency 
between n^, rig, and n^, we have 



/ A -/lR,D „d I ^B,D 



.,, . nj + . b.d ng ~ CM(o, N I nxn ) 

n' 2 = <~CAA(0,iVoInx„) 
Now, if we use the MRC method for detection of x A , we obtain 

r A 12 , I^b,d^r,d| 2 n . / (^r,d^b,d)* / . ,* / n 

yA = (IVo| + , |2 , | 2 ) x a + (t/; i2~TTI ^ n i + ^A, D n 2 ) (15) 

|"b,d| + |"r,d| I"b,d| + I"r,d| 

A similar procedure can be employed for x B by multiplying two sides of equation dTQb by another matrix representing 
the null space of the vector, conveying user A's signal (similar to that of user B). Then, the SNRs at the outputs 
of MRC receivers can be written as 

snRad = (IVdI 2 + \ h [ h ^+]hL\-)w Q p , 16) 

SNR BD = (|fe B , D P + | ^W 2 )iP 

Computing these factors, the destination node chooses to decode the user with the maximum SNR at first. Let u max 
denote the user with the highest SNR, i.e. 

omo rcii, 12 , I^a,d^r,d| 2 a P /,, ,2 . |^b,d^r,d| 2 n P -> 

SNR '<»D = max{(| Vol +77 ^—77 jj) -77", (I | +77 77 r^-rr} (17) 

I"a,d| + I"r,d| iV o I^b.dI + I"-r,d| iV o 

For this user, the destination utilizes the MRC receiver. Then, the effect of this user signal is subtracted from the 
received signal Y RD , i.e. 

Yru - ^M rallx ,DX Umax . 

Then, conditioned on the detection of the user u max being error-free, for the second user (denoted by v? max ) a 
diversity order of two can be achieved, with the following SNR 

SNR„2 D = (\h 2 d | 2 + |Vd| 2 )^ (18) 

max ' ^0 

Corollary III-.l (Generalization). Denoting the set of all users sharing a relay node in a MARC by U, the SNR 
of the user to be decoded ( at the destination ) first is given by 



SNR u i d = max 

max s-n i 



\h I 2 ■ 1 



]KW + W- 

v&A\u 



After decoding the first user's signal, its effect will be subtracted from y RD . Then, the next best user is found, 
similar to equation ([T9T l except the fact that, now we have one user less the initial case. A detailed description of 
the proposed detection algorithm is illustrated in Algorithm [1] 

Algorithm 1 V-BLAST approach for detection of multi-user MARC 

Input: Channel Coefficients (h u p,u £ U and h R , D ), y u o 7 u £ U, yRD 

Output: x„, u £ U 

Initialize: U new <- U, ygg 1 " <- y RD 

for i = 1 to \U\ do 

Find the user with u* nax = argmax 

ue«„. 



|^u,d| 2 h — — i 1 



veu\u 

Decode the signal of u l max , Sc u i using nulling technique for two signals y§§" and y u i D 
end for 

Return x U) u £U 



IV. Performance Analysis in the Error-Propagation-Free state 



In this section, we evaluate the performance of the proposed receiver. We assume that the user-relay channel 
(URC) is stronger than both the user- and relay-destination channels (that is > f2_>D )■ We also assume 

that whenever the relay cooperates, it cooperates correctly. This assumption sounds, since if a strong CRC (Cyclic 
Redundancy Check) code is used, the relay can can detect all users' signals. Hence, if the relay detects an error in 
its received codeword, it discards the whole block and does not cooperate. In the following, we assume a BPSK 
modulation|l For the BEP evaluation, we first compute the probabilities of each state summarized in Table U as 
follows 

po = Pr(s ) = ParPbr 

pi = Pr(si) = (1 - Par)pbr , 2Q) 

p 2 = Pr(s 2 ) = Par(1 - Pbr) 

p 3 = Pr(s 3 ) = (1 - Par)(1 - Pbr), 

where par and (pbr) are the probabilities that the relay cannot correctly decode user A (B) transmitted block. 
These probabilities can be assumed to be negligible at high SNRs. In this section, we derive the BEP for an uncoded 
transmission, and then we generalize our result to a convolutionally coded transmission, in which the user's data 
are first convolutionally coded and then the coded block is transmitted. 

Here we derive the BEP for user it G U in the illustration in Figure [T] Due to the symmetry in the structure, it 
will be the same for user B. We state the main result for the non-coded scheme in the following theorem: 

Theorem IV-.2. Let us denote the BEP for user u £ U in a MARC by Pb, u - Moreover, assume that the relay 
cooperation is error-free. Then, for i.i.d. Rayleigh block fading channels between users (and relay) to destination 
with parameter £1^d, we have 

2- 



Pb.u < 



{P0+P2, 
4 



■D 



■D 



r(2.5)(4fU D +w ) 2 



+ 



F 2, 0.5; 2.5; ^ 



El 
4 

+ I 



O^D \2 



(2.5)(4fU D +w ) 3 



F (3, 1.5; 2.5; I7r ^J — ) 



(21) 

where snr^ : = and F(a,b; c; z) is the generalized Hypergeometric function H26\l and pi's are defined in 
equation 



Proof: Please refer to Appendix |A] ■ 
In the following theorem, we generalize our results for a system which uses the same convolution codes for 
source and relay transmissions. 

Theorem IV-.3. For a Convolutionally coded system using a code with rate r and free distance df ree and k bits 
in each coded transmitted block, when utilizing a hard-input Viterbi decoder, at high SNR, the bit error rate(BER) 
has the following bound 

P b ,u <«c{(|[(pO+P2)jJfc+p/ ''__ ]' 

n 



where a. 



(4fU D +.™f,) 2 T(2.5) 



F(2,0.5;2.5; 



srarS 



^ 4 



+ 



(_U- ± d_ rZ\ 4i2^r> eve 1 c.o c. ™'o 

I^d+wSV' 1 Lr(2.5)(24CJ >D +wS) 3 - r V°' ± ' J ' 4U^ D +.mr 

)"}} 



(22) 



and snr^ : = 



P v rd f 



, , in which ee is the number of codewords with weight df r 



Proof: Please refer to Appendix ICl 



Corollary IV-.4 (Asymptotic behvior). For large SNRs (large values of snr^ or snr^), as snr (snr ) 
have 



lim P bu < lim {a(p +P2) ^— r + (/3pi + 7P3) 

w ->oo w -S>oo SnrQ + A 



1 



snro + A 



+ P30 



where a, (3, and 7 are some positive constants, and g(x) = o(f(x)) means lim ^\ 




00, we 



(23) 



2 The proposed framework can work for any modulation. For simplicity of performance evaluation, we have considered a BPSK modulation. 
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Thus, since at large SNRs, po and pi will be negligible compared to P3, the dominant terms will be the second 
and third terms which decrease inversely with squared SNR. As a result, with the definition of the diversity order, 
the proposed scheme can achieve second order of diversity. 

V. Numerical Experiments 

In this section, we present some numerical results to evaluate the performance of the proposed scheme. We set 
f2_> R = 1 and change f2_>D to consider various link qualities (it will be discussed later). First we investigate the 
uncoded scheme. Figure [3] shows the plots of BEP versus SNR per bit for direct transmission, Alamouti scheme 
using 2 transmitting antennas, and our proposed scheme. For our scheme, the plot based on the analytical result 
of Theorem 1 has been included as well. In this plots, we have assumed ideal user and relay channel (URC). That 
is users to relay links are error free. As we expect, for the proposed scheme, the simulation result coincides with 



10° 
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Fig. 3. The plots of BEP versus SNR per bit for uncoded proposed MARC with ideal URC, uncoded direct transmission, and uncoded 
Alamouti transmission. 

the bound of Theorem 1 at high SNR, and the slope of the plot is the same as that of the Alamouti scheme which 
indicates that the proposed scheme achieves second order of diversity. 

For the rest of the comparisons, we consider coded schemes in which we assume that users and relay use the 
same convolutional code for transmission where each block consists of 50 information bits, i.e., K = 50 . We 
assume that the users A and B apply CRC to their transmitted block in order the relay to be able to detect its 
decoders' errors, also a convolutional code with rate = 1/3 and df ree = 8 which its octal coefficients are [5, 7, 7] 
is considered for the coded transmission. 

First, we show that the upper bound derived in Theorem 2 is tight enough. Figure @] shows the performance of 
our proposed scheme based on both simulation and analytical results in the presence of ideal URC. 

For comparison, the performance of direct transmission (without the relay cooperation) and the Alamouti scheme(for 
two transmit antennas) are provided. We consider equal total transmission powers for comparison. It means that the 
total power of the user and relay transmissions is the same as the power used in direct and Alamouti transmissions. 
As can be realized, the derived upper bound is tight enough. Furthermore, the proposed scheme substantially 
outperforms the direct transmission. Figure |6]shows the plots of BER versus SNR per bit for our proposed MARC 
scheme with different URC's qualities based on both simulations and the upper bound derived in Theorem 2. In 
this plot, for instance the notion "+3dB URC", indicates that the qualities of the sources A and B to the relay 
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Fig. 4. The plots of BER versus SNR per bit for MARC with ideal URC and direct transmission. 
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Plots of BER versus SNR for three-user MARC, two-user MARC with ideal URC, direct transmission, and Alamouti scheme. 
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Fig. 6. The plots of BER versus SNR for different qualities of URC. 

links are 3dB better than the sources to the destination links. From this figure, the upper bound becomes looser 
when the quality of URC decreases. However, it must be noted that the slope of the plot does not change with the 
URC's qualities. 

In Figure [7] 6, we compare the performance of our proposed scheme for different qualities of URC with Alamouti 
scheme and direct transmission (without relay) for the same transmitted powers, based on simulation results. As 
can be generalized, the slope of BER plots for our scheme for different URC's qualities are the same as the 
Alamouti scheme which indicates a diversity of order two, as expected from our theoretical derivations. In our 
further evaluation, we have also considered a MARC system with three users. We have assumed a ideal URC. Fig. 
7 shows the plots of BER versus SNR per bit for MARC with 2 and 3 users. The plots for the Alamouti's scheme 
and no relay cooperation scheme are included for comparisons. As can be realized again, with only one relay, the 
MARC with 3 users achieves second order of diversity. It must be noted that our scheme for 3 users, by 4 antennas 
can achieve second order of diversity for all users in contrast of the traditional schemes in which 6 antennas are 
required for the same number of users to achieve the diversity order of 2. 

VI. Conclusion 

We have introduced a novel scheme for MARC, which utilizes the network coding in an analogue space. Using 
the nulling detection structure applied in V-BLAST, destination can extract the users' signals from two independent 
received paths. We have evaluated the performance of the proposed receiver and provided a tight upper bound for 
its bit error probability. Our analysis for i.i.d. block Rayleigh fading channels have indicated that the proposed 
scheme can achieve a diversity order of two. We have also provided a tight upper bound of bit error rate for the 
coded transmission in which the users and relay use the same Convolutional code. As have shown for n = 3 our 
proposed scheme can be generalized for a MARC with more than two users and only one relay, in which all users 
benefit a diversity order of two, in contrast to the traditional non-cooperated schemes in which 2n antennas are 
required for the same order of diversity, where n is the number of users. 
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Fig. 7. The plots of BER versus SNR for different qualities of URC, direct transmission, and Alamouti scheme. 



Appendix A 
Proof of Theorem IIV-. 21 

In what follows, without loss of generality, we find the BEP for user A. Then, considering the j— th bit of user 
A's codeword, the BEP may be expressed using the conditional bit error probabilities as follows 

3 

Pr(x A ,j / x A ,j) = e[ ^2 Pr (^A,i x A,j\si) Pr(sj) , (24) 
i=0 

in which Sj stands for the rely's state, are defined in TabelHand the expectation is taken on channel gains. As we 
mentioned in previous section, bit error probabilities of the system conditioned on so and s 2 are the same for user 
A, therefore we have 



Pr(x A ,j / x Atj \S = s ) = Pr(x AJ / x AJ \s 2 ) 
Then, for si, as there exists also relay transmission, we have 

Pr(x A ,j / x A j\si) = E 



Q(y/2\hAv\ 2 P/N 



Q(V2(\h KD \ 2 + \h R , D \ 2 )P/N ) 



(25) 



(26) 



For state S3 , depending on whether or not the user with the highest received SNR at the MRC receiver's output 
( u max defined in equation (fTTT ), or defined in equation ([T8T)), is user A, we will have different values for BEP. 
Because of the symmetrical structure of the scheme we have 



Pr^Lx = A) = Pr(u: 



2 

max 



A)= 1 - 
' 2 



As a result, we simply obtain 



Pr(aU / x A \S = 83) = E(-Q(^2SNR <o J + -Q(^2SNR <a J), 



(27) 



(28) 
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where from equations (fTTT i and ([T8T l, we have 

SNR <os = (|Vd| 2 + \h R p\ 2 )P/N 

Now, we compute equations . Plugging Q(x) < \e~~ in equations d25l ), d26l ), and d28l ), and taking the expectation, 
we obtain 

Pr(x A ,j / ^Aj|so) = Pr(x A ,j / ^A,j|s 2 ) < ^0|Vd| 2 ( p / n o) (30a) 
Pr(x Aii ^ x Aii |si) < ^{P/Nq) (30b) 

Pr(£ A)i / x A>j \S = s 3 ) < ±[<l> x (P/N ) + <i>y{P/No)] (30c) 

where <j)\h AD \*> 4>x(-), and (/>y(.) are the characteristic functions of |/i Ai d| 2 , X, and Y, respectively, which X and Y 
are defined as 

x = |/i A ,d| 2 + |^r,d| 2 ,~ n 

A I j. 12 , 1 ^ 

y = IVdI + z , 

wnere z - \ h ^ ? +\ h ^ ? - 

Due to independency between |/i a ,d| 2 and z, we simply have 

MS) = 0\h A , D \*(s)<j) lhRDl 2(s) 

<j)y{S) = (j) lhADl 2 +z (S) = (j) lhADl 2(S)(j) z (S). 

Then, since |/i a ,d| 2 and |/ir,d| 2 are exponentially distributed with parameter Q^.d, their characteristic functions are 
given by 

<Wl 2 ( S ) = Mr»\ 2 ( S ) = s • ( 33 ) 

For deriving the characteristic function of z we need to know its PDF, so we propound the following lemma. 

Lemma A-.5. Suppose that u and v are two independent exponentially -distributed random variables as follows 

u ~ Aexp(-A„u) 
v ~ A exp(— X v v). 

Define z = -^r-. Then, we have 

J u-\-v ' 

f z {z) = 2\e- 2Xz {2\zK 1 (2\z) + 2XzK (2zX)}, z>0 (35) 

where Kq and K\ are modified Bessel functions of second kind. 

Proof: Please refer to Appendix |B] ■ 
We use Lemma IA-.5I with u = |/j a ,d| 2 and v = |/ib,d| 2 - Consider the following equality ll26l 

foo zl i-l e -az K (Bz)d,Z ~ r ^+b)r(^-b) rp( , V b , 1. , 1. q-l-, 

J z e i\ b \pz)az — ^ a+ py +b f(^+|) r ^ + °> ^ 2 > M + 2 ' pg) 
i?e(/i) > |i2e(6)|,i2e(a + /3) > 

where -ftT„(x) is the modified Bessel function of second kindQEmploying equations (|35T)- (|36T ). the characteristic 
function of Z can be directly derived, as follows 

M S) = 4^ D ^^ D ) T(3)r(l) s 

K ' ^ D (4n^ D + s) 3 r(2.5) v ' 'An^D + s' 

Substituting equations (l33l and (l37l) in (l32l ). <p x {S) and 4> y (S) will be derived. Then consequently by substituting 
these results in equations (I30at - (l30cb and by using equation (l28l) . the desired bound is simply obtained. 



l K v (x) is the solution of second-order ordinary differential equation 

dx 2 dx 
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Appendix B 
Proof of Lemma IaTsI 

We write 

F z {z) = Pr(^ <z) = ¥t{uv < z(u + v)) = j Pr(uv < z(u + v)\v = 0M0<% 

Then, because of independency of it and i>, we have 

Pr(u < -—?-\v = £) = Prfit < -^—) = 1 - e~ A "^ 

Pr(u > rr-^-lv = 0= Pr > T - ^ — ) = 1 

z £ - z 

We may then write 

Pr(^<.) = // =0 A,m+ r^ 2 A,e-^(l-e- A "^K 
= (1 - e" A " 2 ) + A e- A ^(l - e~ K f^)di 

Then, after rearrangement, we obtain 

Pr(^ <z) = (1 - e~^) + J£ A„er^ - f™ z X v e-^ + ^di 

= (1 _ e -A„ Z ) + e -A„ 2 _ J- q Awe -(A.^+A.(C+*)) de (38) 

= 1 _ Aue -(A„*+A„,) f- Qe -(^(f )+A„(0)^ 

We use the following equality ||26l 

/°° x^e'i-^dx = 2(£)$K v (2y/fa), [Beifi) > 0, i?e( 7 ) > 0] 

7 

where K v (x) is the modified Bessel function. So by setting 

r/ = = z 2 A n ,7 = A„ 

and by A„ = A u = A, we easily obtain 

Pr(-^ < z \ = i _ 2Aze- 2zA ^i(2zA) 
u+v ~ ' v y 

leading to 

F z (z) = 1 - 2Aze- 2A ^i(2zA),z > 



Appendix C 
Proof of Theorem IIV-. 31 

Bit error rate (BER) in a BSC channel with error probability p where users use a convolutional code at high 
SNR is given by ETl 

Pb~^B dfm 2 d '- p ^ (39) 

Similar to the proof of Theorem IIV-.2I we compute the BER conditioned on the 4 states Sj's of the relay node. 
Briefly, we compute the conditional bit error probability for case 3. By substituting equation (|28T ) in (l39l we have 

(40) 



P b \ S3 = \B dfm 2 rf ^E{ \ (Q( j2Sm~T a J) ^ + i(Q(J2SNR <a J)%"} 



13 



1 —SZ 

Using Q(x) < keT~ we have 



n, S3 < ^B^f-Ei^e-^^ + -i- e - SNR ^^} (41) 

By the same way, the conditional BEP values for cases 1 and 2 can be computed by substituting equations (|25T ) 
and (1261 ) in (l39l).Then by substituting the conditional BER in (l28l) remaining of the proof is like the proof of the 
Theorem 1. 
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